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New type of quantum kinetic equations of the Fermi particles are derived. The Bogolyubov's type 
of dispersion relation, which is valid for the Bose fluid, is disclosed. Model of neutral Rose atoms 
in dense strongly coupled plasmas with attractive interaction is discussed. A set of fluid equations 
describing the quantum plasmas is obtained. Furthermore, the equation of state of a degenerate 
Fermi plasma is derived. 



In recent years a huge number of works have been devoted to the investigation of collective behavior of quantum 
plasmas using a set of hydrodynamic equations. However, to the best of our knowledge, in the literature the derivation 
of the correct quantum fluid equations is missing. It should be emphasized that it is not possible to derive the fluid 
equations with the pressure and the quantum term, which is called the Bohm's or Madelung's potential, simultaneously 
by the Wigner quantum kinetic equation [1]. Strictly speaking the Wigner equation is not suitable for the construction 
of the fluid equations. 

In this Letter, we derive a new type of quantum kinetic equations of the Fermi particles of various species. Our 
consideration is based on general statements of quantum mechanics, especially on the quantum statistics. We start 
our investigation with a single fermi particle and for this purpose, we employ the non-relativistic Pauli equation [2], 
which reads 
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where = ^'Q(r, i, (?) is the wave function of the single particle species a, having the spin s = 1/2(7 (cr = ±1)- 
A(r, and f(f,t) are the vector and scalar potentials, respectively. The last term in Eq.(l) is the potential energy 
of magnetic dipole in the external magnetic field, the magnetic moment of which is 
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where /i^ is the Bohr magneton and a is the operator of the single particle [2, 3]. 
Use of the Madelung representation [4] of the complex function '^^ 
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where aa{r,t^(j) and Sa{r,t,(j) are real, in the Pauli equation (1), yields the following two equations 
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Equation (4) has an obvious physical meaning. Namely, = \^a^ is the probability density of finding the single 
particle at some point in space with a spin s. Whereas, "Pa = VS* — ^A is the momentum operator of the particle. 
Note that if there is no spin dependence of the wave functions (r, t), then pQ(r, t) becomes the ordinary momentum 
of the particle. It should be emphasized that in Eq.(4) the second, quantum Madelung, term describes the diffraction 
pattern of a single electron [5] . 

In the case when there are no external electric and magnetic fields (i? = 0, = 0), after the linearization of Eqs. 
(4) and (5), we get the frequency of quantum oscillations of a free electron 
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Based on the diffraction pattern of electrons, Born has given a statistical interpretation of the wave function, which 
states that in every point in space at a given time the intensity of the de Broglie waves is proportional to the probability 
of observing a particle at that point in space. Therefore, we shall introduce a density of probability distribution fs 
in phase space for the single particle 

1*1' = / d^^p fs{f,p,t) . (7) 

This function fs, obviously must satisfy the normalization condition over all phase space 

J dPr J d?p fs = J <fr = 1 . (8) 

Moreover, Bogolyubov has introduced a one particle distribution function [6] for the system as a whole from Liouville's 
theorem regarding the distribution function f^{f,Ti.T2...T]\[), (where Tq, is the set of coordinates and momentum 
components for the a particle), and derived the Vlasov and Boltzmann equations in the gas approximation, which 
means that the plasma parameter r] (representing the ratio of the average potential energy < ?7 > of particles 
interaction to the average kinetic energy < >) must be less than unity, i.e., 77 = <C 1. Note that the one 

particle distribution function /(r, p, i) is normalized to unity, whereas the Liouville's function to total number of 
particles, i.e., = Nf{r,p,t). The same relation between and f{r,p,t) in an alternative description of kinetic 
theory was obtained by Klimontovitch [7, 8]. To make it more lucid, we shall give a simple explanation about a single 
particle and one particle distribution function [3, 8]. Namely, the probability density of the single particle is one 
particle per unit volume, n^{r,t) = |^(r, f)p, with dimensions 1/V. Whereas the one particle distribution function 
means that in spite of the large number of particles in the unit volume all of them have only one r and p. This permits 
us to express the number density of particles per unit volume as n{r,t) = J d^p f{f,p,t) = N/V. Therefore, we can 
write 

n(f, t) = Nn'if, t) = iV|*(f, t)\^ = J d^p f{f,p, t) . (9) 

Thus, n(r, t) is the density of quantum particles per unit volume. We have assumed that the total number of particles 
of each species is conserved. 

Non-equilibrium states of a Fermi quantum gas are described by the one particle distribution function fa{f,p, t, a), 
which satisfies the quantum Boltzmann equation 

dfa{r,P,t,a) (dr „\ „ ^, , dpa dfa{f,p,t,a) 



m ■ \dt -y ^"v . 

Equation (10) for quasi-particles in a Fermi liquid was written by Landau [9, 10]. Here C{fa) is the collision integral, 
which describes the variation of the distribution function due to particle collisions, the derivative is determined 
by the force acting on the particle, the expression of which is given by the equation (5). When the spin of particles 
is taken into account, the distribution function becomes an operator with respect to the spin variables a. In this 
case the total number density of particles (r , t) equals 

^^;^/a(f,p,t,CT) . (11) 

If there is no spin dependence of the distribution function, then /„ becomes the ordinary quasi-classical distribution 
function fa{r,p,t). Note that the condition for quasi-classical motion is that the de Broglie wavelength H/pp of the 
particle must be very small compared with the characteristic length L, over which f{f,p,t) varies considerably. 

If collisions between particles were entirely negligible, each particle of the system would constitute a closed sub- 
system, i.e., one could neglect the collision integral in Eq.(lO), and the distribution function of particles would obey 
Liouville- Vlasov equation. 

Hereafter, we consider the system as a spinless. We substitute the equation of motion of single particle (5) (neglecting 
the last term) into the kinetic equation (10) taking into account the definition of the density of particles (9) to obtain 

-g^ + iv.V)U + e^[E+ ^— j — + ^ V-^ AV^ ^ = CiU) . (12) 

It should be emphasized that this is a novel equation with the quantum term, which contains all the information on 
the quantum effects. We specifically note also that this equation is rather simple from the mathematically point of 
view. 
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As an example, we employ the above equation (12) to study the propagation of small longitudinal perturbations 
{H = 0, E = — V<p) in an electron-ion coUisionless plasmas. For a weak field, we look for the electron and ion 
distribution functions in the form /„ = /„o + ^/a, where /„o is the stationary isotropic homogeneous distribution 
function unperturbed by the field, and Sfa is the small variation in it due to field. After linearization of Eq.(12) with 

respect to the perturbation, we assume 6 fa and Sep vary like exp{k ■ f— u)t). 
Using the Poisson's equation 

and assuming the Fermi degeneracy temperature Tp = jf^ [Kb is the Boltzmann coefficient, the Fermi distribution 

2 

function is the step function f^o = Q{eFa ^ s), where epa = "'°2^°' ) much more than the Fermi gas temperature, 
then we obtain after some algebra the quantum dispersion equation 



Fafc^tip 1 2kvFa IjJ — kvFa 

where 
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and w can be more or less than kvFe- Note that for u; ^ fewj^i, « 1. 

Let us first consider the electron Langmuir waves, supposing that the ion mass — > oo and lo ^ kvpe, or the 
range of fast waves, when the phase velocity exceeds the Fermi velocity of electrons. In this case we get the dispersion 
relation 



o o 3A!^f p„ Ti. k^ . . 



which has been previously derived by Klimontovich and Silin [11]. The expression (15) exhibits that the high frequency 
oscillations of electrons of a degenerate plasma remain undamped in the absence of particles collisions. Note that the 
Landau damping is also absent, since according to the Fermi distribution there are no particles with velocities greater 
than the Fermi velocity which could contribute to the absorption. 

As k increases, Eq.(15) becomes invalid, but still ui > kVpe and the Landau damping is absent. We now introduce 

the Thomas-Fermi screening wave vector = -^^^^^ In the limit A;^ A;|,y, w tends to to^e (at — > oo) and 
we obtain from (15) 

u = kvFe (l + 2 exp{ J^' }) . (16) 

If we neglect the quantum term ujq in Eq.(16), then we recover waves known as the zero sound, which are the 
continuation of the electron Langmuir wave (15) into the range of short wavelength. Thus the expression (16) 
represents the quantum correction to the zero sound. 

Special and very important case in a quantum plasma is an one-fluid approximation. We further assume that the 
quasi neutrality 

rie = rii (17) 

is satisfied. This equation along with the equation of motion of ions and the equation giving the adiabatic distribution 

of electrons allow us to define the potential field. In such approximation the charge is completely eliminated from the 
equations, and the Thomas-Fermi length tft = J^J disappears with it. 

In order to construct the one-fluid quantum kinetic equation, we neglect the time derivative in the equation (12) 
of electrons, as well as the collision terms, suppose E = — V<^ and H = 0, and write the dynamic equations for the 
quasi- neutral plasma (17) 

(,.V)/. + v(e^+^-^AVH)f =0, (18) 
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^ + (^.V)/,-Ve^^=0. (19) 
In Eq.(19) we have neglected the quantum term as a small one. Note that the Fermi distribution function of electrons 

fe = r \ > (20) 

exp < — ^ ^ + 1 



satisfies the equation (18). Here U = etp + -^^^^l^^Jn , and /Ue is the chemical potential. 

For the strongly degenerate electrons, i.e., Te — > {p,^ = ep), the Fermi distribution function becomes the step 
function 

/e = e(.. + t/-|!-), (21) 

which allows us to define the density of electrons (rig = n, = n = 2 J j^^^jsf) 

We now express ey from the equation (22) and substitute it into the kinetic equation (19) to obtain 

|,,..V,;-v{e..(^)"'-|--LA^}|=0. (.3, 

This is the nonlinear kinetic equation of quantum plasma in the one-fluid approximation, which incorporates the 
potential energy due to degeneracy of the plasma and the Madclung potential. 

To consider the propagation of small perturbations / = fa{p) + Sf{f,p,t) and n ~ uq + Sn{r,t), wc shall linearize 
Eq.(23) with respect to the perturbations, look for a plane wave solution as gj-i^^-r-^t) ^ derive the dispersion 
equation, which resembles the Bogolyubov's dispersion relation in the frequency range kvpe '> ^ kvpi, 




= J'-^ + ^ . (24) 
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We specifically note here that this type of spectrum (24) was derived by Bogolyubov for the elementary excitations 
in a quantum Bose liquid, and created the microscopic theory of the superfluidity of liquid helium [10, 12]. 

How one can explain the similar spectrum (24) in the Fermi gas? In the dense and low temperature plasma the 
formation of bound states is possible due to the attractive character of the Coulomb force [3, 13, 14, 15]. As is well 
known, when the density of particles increases and the temperature goes to zero, the nuclear reaction leads to the 
capture of electrons by nuclei. In such reaction the charge on the ions (nucleus) decreases. Because we have assumed 
the quasi- neutrality, we have therefore supposed that all electrons are in the bound state with ions. Thus one can say 
that in the one-fluid approximation the Fermi plasma may become the Bose system due to the bound state or this 
approximation may imply the formation of the Bose atoms. 

The same Bogolyubov's type of spectrum (24) follows from Eq.(14) in the range of intermediate phase velocities 
kvpe ^ oj ^ kvpi. In this case, except the expression for the real part of frequency (24), we get the imaginary part 
of to from the dispersion relation (14) 

Imcj = -^k^, (25) 

12 TOi ^ ' 

which is much less than Reco. The expression (25) indicates that in the absorption of oscillations the electrons play 
role, since their random velocities greatly exceed the phase velocity. So that the damping rate is determined by the 
electrons alone. 

We shall now derive a set of fluid equations. The Boltzmann and Vlasov type of quantum kinetic equations (12) 
and (23) give a microscopic description of the way in which the state of the plasma varies with time. It is also well 
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known how the kinetic equation can be converted into the usual equations of fluids. Following the standard method, 
we can derive the equations of continuity and motion of macroscopic quantities from Eq.(23) 

— + V-{nu)=0 (26) 

— + {u-V)u = £L^V — 2/3 + - V^A^/^, 27 

at rui no 2niemi y/n 

where u{f, t) is the macroscopic velocity of the plasma 

Obviously from Eqs.(26) and (27) after linearization follows the Bogolyubov's type of dispersion equation (24). 

The more general set of fluid equations for a kind of particles we can obtain from the equation (12) taking into 
account that any elastic scattering should fulfill the general conservation laws of particle number, momentum and 
energy 

^+V-(n„wJ = 0, (29) 

d<Pa> ,^ ^ /ri 1- fi\ 1 w 1 



+ {Ua -W) <Pa >= ea (E + X b) - — VP„ + -^V^La^/tI^ 

V C / ria Iiria 



1 f <rp 



Moreover, the equation of state of a degenerate Fermi plasma can be derived from the kinetic equation (12) by 

2 

multiplying it by integrating over the momentum and employing the equations of continuity (29) and motion 

(30). In the non-relativistic limit the result is 



^ nj^ ~ "a < £a > y (27rn) 



In = Z T / TTT^ C'(/«) , (31) 



where £„ = 2^"^^ and < >= ^ / ^^^^^3 is the internal energy of particles, and the collision terms can 

be positive or negative. 

To summarize, we have obtained a novel kinetic equation for the Fermi quantum plasma. In our formulation of 
kinetic, as well as hydrodynamic equations the quantum term, which contains all the information on the quantum 
effects, is incorporated. It should be noted that the above kinetic equation has mathematical beauty. This equation 
was used to study the propagation of small longitudinal perturbations in an electron-ion collisionless plasmas, deriving 
a quantum dispersion equation. We have discussed the quantum correction to the zero sound. For the special interest 
case we have derived a quantum kinetic equation in the one-fluid approximation and have shown that, when a charge 
and the Thomas-Fermi screening length disappear from the equation, the solution of a linear kinetic equation leads 
to the Bogolyubov's type of dispersion equation, which is valid in the Bose fluid. It is c;lear that in the dense strongly 
coupled plasmas with attractive interaction the formation of the bound state is possible, so that a neutral Bose atoms 
can be created. A general set of fluids equations, describing the quantum plasma was obtained from the new quantum 
kinetic equation. Finally, we have derived a equation of state of a degenerate Fermi plasma in the non-relativistic 
limit. These investigations may play an essential role for the description of complex phenomena that appear in dense 
astrophysical objects, as well as in the next generation intense laser-solid density plasma experiments. 
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